Prediction of Room Temperature High Thermoelectric Performance in n-type 
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First principles calculations are used to investigate the band structure and the transport re- 
lated properties of unfilled and filled Ad skutterudite antimonides. The calculations show that, 
while RhSb3 and p-type La(Rh,Ru)4Sbi2 are unfavorable for thermoelectric application, n-type 
La(Rh,Ru)4Sbi2 is very likely a high figure of merit thermoelectric material in the important tem- 
perature range f50-300 K. 
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The widespread application of thermoelectric (TE) 
coolers has been hindered by the fact that these devices 
have limited efficiencies related to the lack of sufficiently 
high performance materials to use in the active elements. 
The TE performance of a material is characterized by a 
dimensionless figure of merit ZT = aS 2 T/ (kr,-\-k e ) where 
a and k e are the electrical and thermal conductivity, S 
is the Seebeck coefficient and ki, is the lattice thermal 
conductivity. Almost all TE coolers are based on alloys 
of Bi2Te3, which has maximum ZT values slightly ex- 
ceeding 1.0 around room temperature, but substantially 
decreasing performance as temperature is reduced.!!! 

Recently, there has been a renewed effort to find better 
materials in the temperature range 100-300 K with a pri- 
mary focus on novel systems. This effort has been par- 
tially successful in that three new high ZT compounds 
have been discovered: /3-Zn 4 Sb3p and two filled skut- 
terudites CeFe4Sbi2,S and La(Fe,Co)4Sbi 2 .u While these 
materials all have their maximum ZT values substan- 
tially above room temperature and are not applicable to 
cooling applications, considerable understanding, partic- 
ularly of skutterudites, has emerged leading to the hope 
that even higher performance compositions may be found 
and perhaps that lower temperature operation may be 
possible. Here we report first principles results that when 
combined with existing knowledge of skutterudites ma- 
terials, point to a specific composition with a strong po- 
tential for high values of ZT (in excess of 1.0) over a 
wide technologically important temperature range from 
150-300 K. 

The conventional cubic cell of a binary skutterudite, 
for instance CoSb3, contains eight molecular units (space 
group im3). The metal can be also Rh or Ir, and Sb an 
other pnictogen element (P, As). The metal atoms form 
a cubic sub-lattice partially filled by six almost square 
pnictogen rings oriented according to the cubic direc- 
tions. This results in six bonds connecting one Co with 
the six surrounding Sb whereas every pnictogen forms 
two strong covalent bond with other Sb and two par- 
tially ionic bonds with Co. The high values of ZT, with 
a rare earth filling the empty octans of the cubic cell, 



rests on two important effects: a low glass-like kr, due 
to the effect of the filling atom vibrations Ja and on high 
power factors aS 2 related to the electronic features. In 
La(Fe,Co)4Sbi2, in particular, La-filling has the impor^. 
tant effect to decrease kr, about one order of magnitude,!!! 
with respect binary CoSb3 and to lower the top of the 
valence band to near a low-energy high effective mass 
band. This is responsible for the high thermopower in 
p-type materials .0 

In this letter, we investigate by means of first princi- 
ples calculations the electronic properties of skutterudites 
based on RhSb3 and find that La(Ru,Rh)4Sbi2 is promis- 
ing for n-type TE application. Although the thermal 
conductivity of RhSb3 is not glass-like,Ll we can expect 
that filling with La will yield a low kh material because 
the volume available to La's vibration allows consider- 
able amplitudes and a clear relation between filling and 
fcj, in skutterudites is known experimentally. This con- 
jecture is supported by the experimental finding that kr, 
in LaFe4Sbi2 is lower than 20 mWcm _1 K _1 . Based on 
experimental trends of fcz, vs. void size in skutterudites 
we expect kr, for La(Ru,Rh)4Sbi2 to be in the range 10- 
15 mWcm^K- 1 . 

Our electronic structure calculations were performed 
using a well converged (Rkm = 7.0, Rmt = 2.50, 2.25, 
2.25, 2.20 for La, Ru, Rh, Sb respectively) linearized aug- 
mented plane wave (LAPW)|-,methodJ3 including s, p, 
and d local orbital extensions,^! in the framework of local 
density functional theory (Hcdin-Lundqvist parameteri- 
zation) . A fully relativistic calculation was used for core 
states whereas the valence states were done in a scalar rel- 
ativistic scheme. The Brillouin zone (BZ) sampling in the 
self consistent calculation is done by means of a (4,4,4) 
special point grid whereas the density of states, N(e), and 
related quantities are integrated on 295 k points tetrahe- 
dral mesh. 

We used the experimental lattice constant for LaRu4 
Sbi2 (ar, — 9.26 A); the two symmetry unconstrained 
structural parameters related to the Sb position were de- 
termined by means of a best fit to a quadratic curve of 
total energy calculations exploring 2-parameter surface. 
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We obtain x — 0.1583 and y — 0.3412 in the notation of 
Rcf. [h]. These structural parameters were used also for 
the calculation on La(Ru,Rh) 4 Sbi2 performed in virtual 
crystal approximation (VCA), realized in LAPW adding 
one electron in valence and increasing the atomic num- 
ber of the transition metal of one fourth. As discussed 
below, excellent rigid band behavior was found between 
LaRri4Sbi2 and VCA La(Ru,Rh)4Sbi2 supporting the va- 
lidity of VCA, a posteriori, and implying only weak scat- 
tering by alloy disorder. 
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FIG. 1. Scalar relativistic band structure for RhSb3 in an 
energy window centered at the valence band top. Notice the 
zero-gap character due to 3-fold degenerate band at BZ center. 

RhSb3, whose experimental structural parameters are 
cll = 9.23 A and (x,y) — (.1517, .3420), is a zero-gap 
semiconductor because the top of the valence band at BZ 
center belongs to a 3-fold degenerate representation (see 
Fig. |l|) . Anyway the highest occupied band shows strong 
Sb-p component- that might interact, after filling, with 
La-/ resonance,B as happens both in cobalt phosphide 
and antimonide. 

The band structure (BS) of the corresponding filled 
material, LaRu4Sbi2, is drastically different from bi- 
nary RhSb3 because the effect of La is to introduce /- 
resonance in conduction band (CB) and to make avail- 
able extra-electrons that modify the electronic structure. 
This has been proved studying the BS of hypothetical 
YRri4Sbi2. The different arrangement of the bands near 
ep in RhSb3 with respect of CoSb3, as in Ref. HI could 
be strongly related to the different electro-negativity of 
Rh (Pauling electro-negativity: ep — 2.28) and Co (ep = 
1.88) with respect of Sb (ep = 2.05) inverting the ionicity 
of bonds. 

Undoped LaRu4Sbi2 is a metal but its BS is very 
promising for n-type TE application if alloying can be 
done that brings ep to the gap without changing the 
main features of the heavy multi-valley structure above 
the indirect gap in the CB; p-type TE is not favorable 
at reasonable doping because of the lack of heavy hole 
bands near the Fermi level. 
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FIG. 2. Comparison between the scalar relativistic bands 
for LaRu 4 Sbi 2 (gray line) and VCA La(Ru.7 5 Rh. 25)48^2 
(black). The energy zero is fixed at the respective tops of 
the valence band. Due to the rigid behavior of the BS, the 
scattering effect related to the alloying must be quite small. 

In Fig. H we compare the BS of pure LaRu4Sb 12 and 
La(Ru 75RI1 25)4Sb 12 showing that, after alloying, an in- 
direct gap semiconductor (E l g nd — 0.16 eV) is recovered 
without great qualitative changes. The rigid behavior of 
the BS indicates only weak electron scattering due to al- 
loying and the potential for high mobility. The high effec- 
tive mass multi-degenerate minima in conduction band 
(Tab. |) underline the high Seebeck coefficient that we 
have calculated from the BS. 

Starting from,XAPW BS we derived, using standard 
kinetic theory,E3Ej the temperature dependence of the 
thcrmopower at different doping levels giving a direct 
way to compare with experimental measures. The trans- 
port coefficients are calculated using (atomic units): 

a - = ^fde{~ d -^)N{e) v \el (1) 
S= £fJ M-^)N(ey(e)(s - ep), (2) 

where ep is the chemical potential, e the electron charge, 
fo the Fermi distribution function and v(e) the average 
velocity of electron with energy s; the constant scattering 
time (r) approximation is assumed. 

The measured values of Seebeck coefficient for 
LaFe 3 CoSbi2 and CeljsaSbia are about 100 fiVR- 1 
at room temperaturepfj in the compound we present, 
La(Ru. 75 Rh 25) 4 Sbi2, S ranges from -150 to -200 [iVK - 1 
for different doping from 4.5 10 18 to 3.7 10 20 electron 
cm~ 3 . As shown in Fig. |^ the temperature dependence 
is very favorable having a maximum below 200 K and 
remaining quite stable at higher temperature. 
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FIG. 3. Temperature dependence of Seebeck coefficient 
for La(Ru.75Rh.25)4Sbi2 at different doping concentrations in 
electron per cm -3 . 

In summary, we have presented band structure calcu- 
lations for RhSb3, LaRu4Sbi2 and La(Ru.75Rh. 25)4^12. 
Our results and the experimental trends suggest 
La(Rui_ K Rh :E )4Sbi2 as a very promising material for n- 
type TE application because of (1) expected low thermal 
conductivity, (2) multi- valley high effective mass char- 
acter of the CB, (3) rigid behavior of the band struc- 
ture after alloying, (4) high Seebeck coefficients in the 
temperature range useful to cooler device applications. 
It should be emphasized that although the results indi- 
cate weak n-type carrier scattering due to Ru-Rh dis- 
order, the sensitivity of the band structure to La filling 
indicates strong scattering by La defects, so preparing 
samples with maximum filling is key to assessing the TE 
potential of La(Rui_ a: Rh a: )4Sbi2. We underline that the 
difficulty to prepare n-type highly filled skutterudites has 
been overcame using high pressure synthesis,li-3 for CoSb 3 
and this should work also here. 

We are grateful to Dr. F.J. Di Salvo for the careful 
reading and comments on our manuscript. This work 
was supported by ONR and DARPA. 
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TABLE I. Effective masses for VCA La(Ru. 7 sRh.25)4Sbi2 
along the high symmetry directions in BZ at F, N and at the 
minimum near A direction, K„ = (0.170, 0.486, 0.001) in unit 
of 2-7r/a. e denotes the band energy with respect to the top 
of the valence band (see FigjiJ); r\ is the number of equivalent 
valleys. The effective masses are along A, A, £ for r"; along 
£, G and D for N c . For the effective masses are along 
the cubic directions. Both the six N c and 24 pockets are 
active in transport above 150 K for any n-type doping. 
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